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Abstract 

Several conditions are given when a packing of equal disks in a 

torus is locally maximally dense, where the torus is defined as the 

quotient of the plane by a two-dimensional lattice. Conjectures are 

presented that claim that the density of any strictly jammed packings, 

whose graph does not consist of all triangles and the torus lattice is the 

standard triangular lattice, is at most -tt-5=, where n is the number 
° ' n+l sj\2 ' 

of packing disks. Several classes of collectively jammed packings are 
presented where the conjecture holds. 
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1 Introduction and definitions 



Packings of disks are of great importance in number theory, granular mate- 
rials, algebraic algebraic number theory, and who knows what else. Their 
rigidity properties are also important and can provide useful tools for their 
analysis, in particular for finding maximally dense packings, or locally max- 
imally dense packings, in certain situations. Here we primarily investigate 
the case of periodic packings of equal disks, where the period lattice is the 
ubiquitous triangle lattice given by the edge-to-edge tiling of equilateral tri- 
angles. One reason for choosing this lattice is that it is the basis for the 
most dense packing of equal disks in the plane. See, for example, the classic 
books [9j H2], for a good discussion of these problems. Another reason for 
this choice is that the techniques of rigidity theory can be applied here, while 
in many other cases some of the tools described here are not available. 

1.1 Lattices 

A lattice A in Euclidean space K d is the set of all integral linear combinations 
of the vectors g 1: . . . ,g D in E d , 

A(fli, ...,g D ) = {nwi H h n D g n | m, . . . , n D <E Z}. (1) 

Note that D can be greater or less than d, but the usual definition of a lattice 
is when D = d and gi, . . . ,go is a basis of E d . This sort of lattice can be 
regarded as a projection, or injection, of the usual lattice into E d . 

We also find it is useful, for d = 2, to use complex notation, so C = E 2 , 
and define g{9) = e ld = cos 9 + ism 9. For notational fluidity, we define 

1 \/3 

9a = #(27r/6) = - + — % and A A = A(l,</ A ), 
go = 0(27r/4) = i and A n = A(l,g D ), 

which is the usual triangular lattice and the square lattice as in Figure [l] 
We say that any lattice of the form A(z, g&z) is a triangular lattice. In other 
words there are generators of equal length such that the angle between them 
is 2n/6. Similarly a square lattice is of the form A(z,g\jz). 

For any two-dimensional lattice A in the plane, not contained in a line, 
we define the corresponding torus as 

T 2 (A) = C/A, 
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where two points are identified if they differ by a vector in the lattice A. 
We regard T 2 (A) as a compact metric space, where locally it is isometric to 
the flat Euclidean plane. Alternatively, we can regard T 2 (A) as obtained by 
identifying opposite edges of a fundamental parallelogram as in Figure [T] 

Since T 2 (A) is compact, it has a finite area, namely the area of any 
fundamental parallelogram. For example, 

Area(T 2 (A($ri,0 2 )) = |#i||# 2 | sin0 = Im^^), (2) 

where 6 is the angle between g t and g 2 , z is the complex conjugate of z, and 
Im(z) is the imaginary part of z. 




Figure 1: This shows a portion of a triangular lattice A a and the square lat- 
tice An, along with the corresponding circle lattice packings. A fundamental 
region is shaded. The circle lattice packing for T 2 (Aa) has the densest pack- 
ing density of any packing in the plane. 

We say that a triangular torus is any torus of the form T 2 (A(z, g&z)). 
In other words, a triangular torus is where the lattice used to define it is a 
triangular lattice. Similarly a square torus is defined by a square lattice. 

1.2 Packings 

A collection of n equal (radii) disks with disjoint interiors in T 2 (A) is called 
a packing. When n = 1, the packing is called a lattice packing. The density 
of a packing is 5, the ratio of the sum of the areas of the disks divided by the 
area of the torus. So the density of n disks of radius r in T 2 (A(z, g(9)z)) is 
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5 = nnr 2 /(\z\ 2 sin 6). 
For the triangular lattice packing when r = ] 



1/2 and n = 1, 



(3) 



5A = 7r(i)7sin(27r/6) 



(4) 



which is the maximum density for any packing of equal disks in the plane. 

When the packing is lifted to a packing of an infinite number of packing 
disks in the plane, this definition of packing density agrees with the usual one, 
where larger and larger regions are intersected with the packing to calculate 
the packing density. 

One of our main problems is to understand the most dense packings of n 
equal disks in a triangular torus. For certain values of n, this determination 
is easy, namely when there is a triangular lattice packing inside the triangular 
torus. This happens when Aa is a sublattice of A(z, g^z). So z = ri\ + n 2 g^, 
and 

\z\ 2 = zz = (n x + n 2 gA){ni+ n 2 g&) = n\ + nin 2 {g& + g&) + n\ 



The ratio of the areas of the tori is given by (J5J), since the angle between the 
generators 2%/6 is the same in both cases. Furthermore, ^ is the number of 
packing disks by ^ since the density is the same for both packings. Putting 
this together, we get the following: 

Theorem 1 The maximum density of a packing of n equal disks 5a in a 
triangular torus is achieved when and only when n = n\ + n 1 ra 2 + n\ for 
rii, n 2 integers. 

Figure [2] shows this for n± = 2, and n 2 = 1 and n = 7. 

We call any integer of the form ^ a triangle lattice number, which are 
also called Loeschian numbers in [5], page 111, Sloane sequence A003136. 
The first few numbers of this form are 0, 1, 3, 4, 7, 9, 12, 13, 16, 19, 21, 
25, 27, 28, 31, 36, 37, 39, 43, 48, 49. It is easy to see that the product of 
triangle lattice numbers is a lattice triangle number, and it is shown in [T8] 
that n > is a triangle lattice number if and only if the non-square prime 
factors are congruent to 1 modulo 6, or it is the prime 3. 



n\ + n x n 2 + n\. 



(5) 
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Figure 2: This shows a triangular packing of 7 equal disks in a triangular 
torus of maximal density n / \fY2. The fundamental region of the sublattice 
is indicated. 

1.3 The conjectures 

One of our motivations is the following conjecture of L. Fejes Toth in [TUtlTT]. 
A packing of circles in the plane is said to be solid if no finite subset of the 
circles can be rearranged so as to form, along with the rest of the circles, a 
packing not congruent to the original one. 

Conjecture 1 (L. Fejes Toth) The triangle packing in the plane, minus one 
packing disk, is solid. 

The definition of a solid packing refers to unlabeled packings, so permut- 
ing the packing elements is not considered a distinct packing. The following 
is a related conjecture, but in terms of finite packings in the triangular torus 
is what we call the Density Gap Conjecture. 

Conjecture 2 Let n a positive integer so that n + 1 is a triangular lattice 
number but n is not a triangular lattice number. Then the most dense packing 
of n equal disks in the triangular torus is 5 = ^t[0~a = ^pi^j; an d is 
achieved only by the triangular packing with one disk removed. 

The point is that Conjecture [2] implies Conjecture [T] since a finite sub- 
packing of a counterexample to the Conjecture [T] can be placed in a triangular 
torus with a sufficiently large number of disks. We do not provide a proof of 
Conjecture |2j but we do provide some evidence for it in the following sections. 

Following the definitions in [5], we say that a packing is collectively 
jammed if the only continuous motion of the packing disks, fixing the lattice 
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defining the torus and fixing the packing radius, is a continuous translation. 
The following is another conjecture closely related to Conjecture [2j 

Conjecture 3 If a packing of n equal disks in a triangular torus is collec- 
tively jammed and is not a triangular packing, then its density is at most 
^pi^A = ^+1^2 ' an< ^ ^ e max ^ ma ^ density is achieved only by the triangular 
packing with one disk removed. 

Note that Conjecture [3] concerns any collectively jammed packing of disks, 
whether or not there is a triangular packing of n + 1 disks available. 



1.4 Results for small numbers of packing disks 

In [7] the locally maximally dense packings of n equal disks in a triangular 
torus were investigated, and all the locally maximally dense packings were 
found for 1 < n < 6. Figure [3] shows those packings. For each of those 
values of n, there is only one locally maximally dense packing, and thus it is 
the maximally dense packing. For all of these packings, except n = 5 these 



n=2 











- 











n=4 



n=5 




Figure 3: These are the locally maximally dense and maximally dense pack- 
ings of equal disks in a triangular torus for n from 1 to 6. A fundamental 
region is shaded in each case. 

are the packings expected from Conjecture [2j namely, either a triangular 
lattice packing or such a packing minus one packing disk. For five disks, a 
fundamental region is shown that has three-fold dihedral symmetry, and thus 
the whole packing has that symmetry. Similar results were found in [B] for a 
square torus. 
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1.5 Rigidity theory 



For any packing the graph of the packing is obtained by connecting each 
pair of centers of the packing disks that are tangent to each other with 
an edge. Indeed the (labeled) disk centers form a configuration of points 
P — (Pi;P2, • • • ) in the plane or quotient torus, and we consider each edge 
{j, k} as a strut in the sense of a tensegrity framework as in (3J, HI [2]. An 
infinitesimal flex of framework G(p) is a sequence of vectors p' = (p' 1; p' 2 , . . . ) 
where each p^- is associated to pj such that for each strut {j, k} 



Although a periodic configuration p is infinite, for a fixed lattice A, we con- 
sider each equivalence class Pj + A as a single point in the quotient torus, 
and most of our calculations are done in the Euclidean plane. The only con- 
tinuous isometries of the quotient torus are translations, so the only trivial 
infinitesimal flexes of our torus tensegrities are when each p ! - are the same 
constant. With this in mind, we say that a strut tensegrity (or more gen- 
erally any tensegrity) is infinitesimally rigid if the only infinitesimal flexes 
are trivial. A strut tensegrity, in a torus, is rigid if the only continuous mo- 
tions of the vertices not decreasing the lengths of edges are translations. The 
following is a useful result that can be found in [TJ [2] : 

Theorem 2 Any strut tensegrity in a fixed torus is rigid if and only if it is 
infinitesimally rigid. 

That infinitesimal rigidity implies rigidity is a general fact. The converse 
implication is because all the edges are struts. 

If the infinitesimal flexes p' that satisfy ^ are equalities, we say that 
p' is an infinitesimal flex of the bar framework. Again if the only infinites- 
imal flexes of the bar framework are trivial, we say the bar framework is 
infinitesimally rigid. 

Any assignment of a scalar Ujk = ojkj to each edge {j, k} of the packing 
graph is called a stress (where Ujk = when {j, k} is not an edge in the 
graph) and it is called an equilibrium stress or equivalently a self-stress, if 
for each vertex p., 



The following result of B. Roth and W. Whiteley [T9] is very useful in 
deciding when a strut tensegrity is infinitesimally rigid, because often solving 




(6) 




k 
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the equality constraints of the bar conditions are somewhat easier than the 
direct inequality constraints of 

Theorem 3 A strut tensegrity is infinitesimally rigid if and only if the un- 
derlying bar framework is infinitesimally rigid and there is an equilibrium 
stress that is negative on all the struts. 

So the lack of an equilibrium stress, all of the same sign, or an infinitesimal 
motion of the bar framework, is enough to show that the disks can be moved 
keeping the packing property. 

With these results in mind we can model the local behavior of a packing 
of disks with a strut tensegrity on the edges connecting the centers of tangent 
packing disks. If the underlying graph of a packing is rigid (and therefore 
infinitesimally rigid) as a strut graph in a fixed torus, then the packing is 
collectively jammed as defined in Subsection |1.3| 

A consequence of Theorem [2] and Theorem [3] is that if a packing of disks 
is collectively jammed, then m > 2n — 2 + 1 = 2n — 1, where n is the number 
of packing disks, and m is the number of contacts, i.e. the number of edges 
in the contact graph. 

We say that a packing is locally maximally dense if there is an epsilon 
e > such that any e perturbation of the packing to another packing (with 
possibly different common radius) has the same density. Another conse- 
quence of the proof of Theorem [2] is that any locally maximally dense pack- 
ing must have a subpacking that is rigid and thus infinitesimally rigid. We 
call such a rigid subset of the packing disks a (rigid) spine of the original 
packing. In other words, the spine is collectively jammed. In many cases the 
spine consists of the whole packing, or there are isolated packing disks not in 
the spine. A packing disk, not in the spine is called a rattler. For example, in 
[8] there are examples of locally maximally dense circle packings in a square 
torus, or a nearly square torus, that has rattlers. 

The following is our main result, which is a special case of Conjecture [2j 
The next sections will provide a proof. 

Theorem 4 If the graph of a collectively jammed packing of n disks in a 
triangular torus consists of only (equilateral) triangles and rhombi, the degree 
of each vertex, which is adjacent to a triangle, is at least 5, and the degree of 
each vertex adjacent to all rhombi is four, then it consists of a single strip of 
rhombi, all rhombi, or all triangles. In these collectively jammed cases, its 
density is 5 < ^6 A = ^j^=. 
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2 Triangle and square tilings 

In this section we provide information about the rigidity and flexibility prop- 
erties of triangle and square tilings of the plane. An example of a square and 
equilateral triangle tiling is shown in Figure |4j 



Figure 4: This shows a portion of a square and triangle tiling of the plane. 

A circle packing is obtained by centering circles of radius one-half of the 
edge length at each of the vertices. In fact, any graph embedding whose edge 
lengths are the same, and such that the internal angles are all greater than 
27r/6 will correspond to a circle packing. 

Let g = e 27 ™/ 12 , a twelfth root of unity in the plane C. g is one of the 
four roots of the cyclotomic polynomial z 4 — z 2 + 1, which can be regarded 
as a defining relation for g, and all 12 powers of g are an integral linear 
combination of l,g,g 2 ,g 3 . So A = A(l, g, g 2 , g 3 ) is a four-dimension lattice, 
and any point z in A is a unique integral linear combination n\ + n 2 g + 
n^g 2 + n^g 3 = z, since the cyclotomic polynomial is irreducible over the 
rational field. (Alternatively, this follows since 1 and a/3/2 are independent 
over the rationals, and they form the real and imaginary components of the 
four powers of g.) 

Position the origin at one of the vertices of the tiling and align one edge 
of graph with one of the twelve generators of A. One sees that each oriented 
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edge in the tiling is one of the powers of g and each vertex of the tiling is in 
A. Also g 2 is a sixth root of unity with cyclotomic polynomial z 2 — z + 1, 
so Ai = A(l,g 2 ) and A 2 = A(g,g 3 ) are both sublattices of A. Each directed 
edge (directed in either direction) of each triangle in the triangle and square 
tiling belong to exactly one of Ai or A 2 . We say that two polygons in the 
plane, with disjoint interiors, are edge- adjacent or edge connected if they 
have an edge in common. We say they are vertex- adjacent if they have 
a vertex in common, but are not edge-adjacent. When two triangles are 
edge-adjacent, their edges belong to the same sublattice, and when they are 
vertex- adjacent to each other and edge- adjacent to the same square, their 
edges are in opposite lattices. 

Consider an infinite edge-to-edge tiling of the plane by triangles and 
squares such that there are no infinite strips of squares or triangles, as in 
Figure |5j 




Figure 5: This shows a strip in a square-and-triangle tiling. 



Theorem 5 The bar graph consisting of the edges of a square and triangle 
tiling of the plane, with no infinite strips, flexes with one degree of freedom 
by rotating the edges of the lattice A% relative to the lattice A 2 . If the tiling is 
periodic with respect to a triangular lattice, it remains periodic with respect 
to a contracted triangular lattice. If the tiling is periodic with respect to a 
square lattice it remains periodic, and the period lattice generators remain 
perpendicular, but they do not remain the same length. 

Proof. If an edge of the square and triangle tiling is in Ai, it must be 
an even power of g, and if it is in A 2 , it is an odd power of g. This is 
indicated by the shading in Figure [6] and Figure [7j Note that both Ai and 
A 2 are two-dimensional discrete lattices, so they can be rotated by gg = e tB 
and the same integral linear combination of the generators of A can be used 
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for Ai and A 2 . If there are no infinite strips, there must be edges of the 
tiling that are part of the generators from both Ai and A 2 . So, explicitly, if 
z = rt\ + n 2 g + n^g 2 + n^g 3 is in A, then the flexed z is 

z(9) = m + n 2 gg g + n 3 g 2 + n 4 g 3 g e = X 1 + X 2 g e , (7) 

where Aj = % + /i3<7 2 is in A 1; and A 2 = n 2 g + n^g 3 is in A 2 . 

The next question is to determine what the lattice of the flexed configura- 
tion is. Let Ai + A2 represent one of the generators of a period parallelogram 
for a triangular lattice tiling of triangles and square, where Ai G Ai and 
A 2 G A 2 . Then g 2 (X\ + A 2 ), the rotated vector by 27r/6, is another generator. 
Then in the flexed configuration these generators are flexed to Ai + X 2 ge, and 
g 2 (Xi + X 2 ge), respectively, since g 2 X± G Ai and g 2 A 2 G A 2 , and clearly the 
length of the these generators are |Ai + X 2 gg\ = |fl ,2 (Ai + X 2 ge)\. But for the 
same pair of generators g 4 (Ai + A 2 ) is the difference vector between them, and 
by a similar calculation it also has the same length as the other generators. 
Thus the lattice is a contracted triangular lattice, as desired. 

For the square lattice, there is an important difference in how the gen- 
erators are flexed. If Ai + A 2 is a generator as before, then g 3 (Ai + A 2 ) is 
the other generator of the same length at right angles. But now g 3 Ai G A 2 , 
and g 3 X 2 G A 1; so the flexed generators are g 3 {Xige + A 2 ) and g 3 (Ai + X 2 go), 
respectively. But nevertheless, as long as there are no strips, the infinite 
framework can still be flexed as before, and when gg = g 2n /e all the squares 
have been flexed to tt/6 rhombi, and the vertices are part of a triangular 
lattice. If the period parallelogram remains a square, it means that square 
lattice is a sublattice of the triangular lattice. But that is not possible, since 
if A G Ai, then its rotation by 27r/4, given by g 3 , would be in A 2 . But these 
two lattices are disjoint except at 0. □ 

Note that a consequence of Theorem [5] is that that any packing whose 
graph consists of triangles and squares, with no strips in a triangular torus, 
is not rigid, i.e. not collectively jammed. 

3 Triangle and rhombus tilings 

We next generalize Theorem [5] to the case when the quadrilaterals are not 
necessarily squares, but rhombi. For the case when the tiling is the graph 
of a packing of equal circles, we assume that each internal angle of each 
rhombus is strictly between 2n/6 = 7r/3 and |27r = 2ti/3. If the angle is 
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Figure 6: This shows a periodic triangle and square tiling with a period 
rhombus shown and the rigid groups of triangles shown. The triangles with 
edges in the same lattice have the same shade. The 3- fold covering of this 
tiling is the one given in Figure [4] 

either of those extreme angles, we assume that the rhombus is subdivided 
into two equilateral triangles. If the graph has that property, we say that it 
is completed. 

Theorem 6 Suppose a (edge-to- edge) tiling of the plane into (equilateral) 
triangles, rhombi, such that each vertex that is adjacent to a triangle has 
degree 5 or 6, and there are no infinite strips. Then all the rhombi are 
congruent and the tiling can be flexed to a square and triangle tiling as in 
Theorem [3| In particular, if the tiling is periodic with respect to a triangular 
lattice, it remains periodic with respect to a triangular lattice. 

Proof. The edge-connected triangle components form convex polygons 
with between 3 and 6 vertices, with internal angles other than tt, and such 
that the internal angle at each other vertex is either tt/3 or tt/6. There are 
no non-convex vertices because of the completed property. Figure [8] shows 
examples of these polygons for all the cases from 3 to 6. 

On each edge of each of these polygons, there is a rhombus or adjacent 
rhombi such that they form a parallelogram edge-adjacent to one edge of the 
triangle component. The external angles of the triangle component are such 
that when you count the angles, not coming from the edge-adjacent rhombi, 
they add to 2tt. In cases (b), (c), and (d) there are 6 such angles, so the 
average angle is 2tt/Q = tt/3. But since each angle must not be less than 
7r/3, they are all tt/3. In case (a) there are three such angles, the average 
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Figure 7: This shows the flexed framework of Figure [6] where each square 
is flexed until it is the union of two triangles. When the extra edge is filled 
in, the resulting tiling is a triangle lattice tiling. There are 13 vertices in 
the original tiling. When it is flexed even further in Figure (c) to collapse 
the squares entirely, another tiling of triangles occurs, in this case with 7 
vertices. 

is 27r/3 and each angle cannot be greater than 27r/3. So all the angles are 
27r/3. Thus each triangle component is vertex- adjacent to another triangle 
component at each vertex (not including the vertices whose interior and 
exterior angles are tt). Furthermore, as one proceeds around the boundary 
of each triangle component, the rhombi are all congruent, with the same 
internal angles (when they are not squares), appearing in the same order 
with respect to a clockwise order, say, along the boundary. So the rhombus 
components are all parallelograms with those same angles, since there are 
no infinite components. So each triangle component has its edge directions 
in one of two possible sets of directions, where each set consists of three 
directions 27r/3 apart. Thus one set can rotate relative to the other. The 
whole graph can be flexed to a triangle and square configuration, and the 
conclusion of Theorem [5] holds. Figure ^h) shows an example of a case when 
the rhombi are not squares. □ 
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Figure 8: The shaded polygons are examples of components of triangles. 
The edge- adjacent rhombi are shown as well. The extra edges for cases (c) 
and (d), when the internal angles are ir/6, are there because of the degree 5 
hypothesis. 



4 Rigid rhombus and triangle strip tilings 

When all the faces of packing graph are rhombi, we provide an explicit de- 
scription what they look like in Subsection |4.1 Then in Subsection 4J3 we 



calculate the maximum density of such packings for a triangular torus. 



4.1 A description of a collectively jammed rhombus 
packings 

Another case of tilings corresponding to strictly jammed packing graphs in 
a triangular torus is when all of the faces are rhombi. We first determine the 
abstract graph of such a packing. For any edge-to-edge tiling of a torus let 
Vj be the number of vertices of degree j for j = 3,4,..., and let fj be the 
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number of faces with j sides, for j = 3, 4, . . . . Define the averages 

1 :\ ^ :\ 

3 J 3 

where t> = t>j is the total number of vertices, and / = fj is the total 
number of faces. (Notice that all the sums are finite since there are only a 
finite number of vertices and faces. The bar notation, used here for averages, 
is not to be confused with complex conjugation.) Let e be the total number 
of edges in the graph. Then counting the edge- vertex adjacencies and the 
edge-face adjacencies we get 

vv = 2e = ff. 

Assuming that each face is simply connected, and since the Euler character- 
istic of a torus is 0, we get 

= v - e + f = 2e/v - e + 2e/f, and 
1/v + l/f = 1/2. 

So when all the faces are quadrilaterals, / = 4, and thus v = 4 as well. But 
if there are any vertices of degree 3 and the graph is a packing graph, the 
three angles at any such vertex must each be 2ir/3, implying that there is a 
triangle in the graph. Thus t>3 = 0, and so all the vertices must be of degree 
4. Summing up, we have shown: 

Proposition 7 Any packing graph of equal circles in a torus, where each 
face is a rhombus, has each vertex of degree four. 

So far we have not assumed that the torus is a triangular torus, or even that 
the packing is jammed in any sense. We next consider when the packing is 
collectively jammed. 

If the graph of a collectively jammed packing consists of only rhombi, 
then in the lift to the plane, there are infinite strips of rhombi, where two 
opposite sides of each rhombus are in a fixed direction. The opposite sides 
of each rhombus are adjacent to the next rhombus in the strip. The other 
pair of sides are each adjacent to another rhombus in another strip of all 
different rhombi, where each rhombus in the second strip is edge-adjacent to 
a rhombus in the first strip. Figure [9] shows an example of such strips. 
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Figure 9: This shows a portion of an infinite rhombus tiling in the plane, 
where an attempt is started to find a reciprocal graph is indicated in dashed 
lines, but cannot be continued. 

By Theorem [2] and Theorem [3j a necessary condition on the graph of a 
packing to be collectively jammed is that there is an equilibrium stress on 
the edges that is strictly negative on each edge. One way of picturing this 
stress is to draw the reciprocal graph, which is another graph whose vertices 
correspond to the faces of the original graph, whose faces correspond to the 
vertices of the original graph, and whose edges correspond to the edges of the 
original graph and are perpendicular to the original edges. This reciprocal 
graph is indicated in Figure[9]with dashed lines. In the case when all the edges 
of the original graph are the same length, the stress on an edge corresponds 
to (minus) the length of the edge of the reciprocal. See [TJ [2] for a discussion 
of this with regard to packings. 

Since the packing is periodic, the stress and the reciprocal extends to the 
infinite packing in the lift in the plane. But each reciprocal polygon in the 
lift must be a convex polygon, which in our case is another quadrilateral, 
and the straight lines in the reciprocal, corresponding to the parallel sides of 
each adjacent strip, must not cross. So this means that the lines are parallel. 
Doing this for all non-crossing strips gives the following: 

Proposition 8 If the graph of any collectively jammed circle packing con- 
sists of only rhombi, then the rhombi are all congruent and translates of a 
single rhombus. 
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Proposition [8] supplies the condition for the graph of the packing to have 
a proper stress that is necessary according to Theorem |3j but we still need 
conditions for the graph to be infinitesimally rigid as a bar graph. 

We now know that each edge of the rhombus tiling is in one of two 
directions. As one proceeds from one vertex to the next on the torus in one 
of those directions the vertices are visited one after the other until one gets 
back to the starting vertex. Suppose that Pi,P2, ■ ■ ■ ,Pk is this tour of vertices 
and it does not include all the vertices on the torus. Then there is a non-zero 
infinitesimal flex p\ = p' 2 = • • • = p' k of the bar graph that is on all the other 
vertices not on this tour, where each p'- is perpendicular to the direction of 
the other edge of the original rhombus. Figure 10 shows an example of this. 
It is easy to check that this is the only infinitesimal flex possible. Hence we 




Figure 10: This shows an example of a periodic flex when the tour of centers 
is not all the vertices of the graph of a rhombus packing. 

get the following: 

Theorem 9 If the packing graph of a circle packing consists of only rhombi, 
the packing is collectively jammed if and only if all the rhombi are congruent 
with the sides in only two directions and the tour of vertices in each direction 
consists of all the vertices in the packing. 

Figure 11 shows an example of the case when the packing graph consists 
of rhombi and each tour in each direction visits each vertex. 



4.2 Strips of rhombi and triangles 

Another case of collectively jammed packings in a triangular torus is when 
there is a single strip of rhombi, and a corresponding strip of triangles of 
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Figure 11: This shows an example of a collectively jammed packing on a 
triangular torus that, where one tour of vertices is indicated. The other tour 
is 1,4, 7,2,5,8,3,6. 



arbitrary thickness in the graph of the packing. This includes, as a special 
case, the situation in Subsection 4A of all rhombi. We call a strip of rhombi 
in a tiling of a torus strip if each the rhombus is edge-adjacent with exactly 
two other rhombi at all the lines containing the intersections between rhombi 
are parallel when lifted to the plane. So these rhombi form a cyclic chain in 
the torus. In the case here, except for case when all the faces are rhombi, 
we assume that there is only such strip of rhombi, and all the other faces 



of the graph are (equilateral) triangles. This is indicated in Figure 12 In 



Section 4.1, we assumed that all the faces were rhombi, and there were only 
two strips of rhombi, each one corresponding one of the two parallel classes 
of edges in the lift. 

Place one of the lattice vectors, the origin, at lower left corner of one of 
the rhombi. Then as you proceed in the direction of the strip, you must come 
back to the same point on the torus for the first time, and thus you must hit 
another lattice point in the lift to the plane. Let a be the number of layers 
of triangles, and b the number of steps needed before you hit another lattice 



point on the upper layer, as indicated in Figure 12 Thus the total number 



of vertices on the torus is n = c(a + 1). The case when a = is the situation 
when all the faces are rhombi as in Subsection 14.11 

In order to do calculations of density, we orient the packing so that the 
strips all are aligned along the real axis, the edges all have unit length, the 
rhombus edge direction, not along the real axis, is given by the unit length 
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c 



Figure 12: This shows a strip tiling of triangles and rhombi, where c is the 
length of the strip of rhombi, a is the width of strip of triangles, and b 
represents the displacement of the identification of the bottom to the top. 
Three of the lattice points are indicated. 

complex number g, where lm(g) > 0, and the internal angle at the origin for 
the rhombus is between tt/3 and 27r/3. The edge directions of the triangles 
are g& and g\ as before. 

4.3 Calculating the density of triangle-rhombus strip 
packings 

Theorem 10 The density 5 of an equal circle packing whose graph is a 
triangle-rhombus strip tiling in a triangular torus is always such that ^pr^A > 
5, where n is the number of packing disks. 

Proof. Since the area of the unit length rhombus is lm(g), the total area of 
the torus is (Im(g) + a\/3/2)c and the density of the corresponding packing 
is 

§ _ n(l/2) 2 n _ 7T c(o + l) _ ?r o + l 
" (Im(#) + aV3/2)c ~ 4 (l m ( g ) + aV3/2)c ~ ±1m{g) + ay/3/2' 

The next task is to use the condition that the packing is in a triangular torus. 
Let v be a generating vector for a lattice that defines the triangular torus, so 
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that lattice is A(v, g&y). Since A(c, g + ag A + b) is the same lattice, we have 

c = n x v + n 2 g A v = (ni + n 2 g A )v 
g + a#A + b = n 3 v + n A g A v = (n 3 + n 4 g A )v, 

where ni,n 2 ,n 3 ,n4 are integers. Then paying attention to the order of the 
lattice generators 



det 



ni n 2 



riiHi — n 2 n 3 = 1. 



Ji3 n A/ 

Then we can solve for v and simplify, 

_ c _ c(rai + n 2 g A 1 ) 
nx + n 2 g A n\ + ri\n 2 + n| ' 

Similarly, we solve for g, 

c(m + n 2 g^){nz + ™40a) 



(9) 



p = -a# A - b + 
= -ag A - & + 



+ %n 2 + n\ 



(10) 



2 Cni^s + ™2«4 + nin 4 5-A + n 2 n 3 #A J. 



nf + ri\n 2 + n 2i 

We are primarily interested in the imaginary part of g, which is 



lm(g) 



—a + 



-a + 



n\ + riin 2 + n 2 



+ riin 2 + n| 



(nin 4 - n 2 n 3 ) 
2 



V3 
2 



by g. Thus, 



Im(g) + a 



V3 



2 nf + rtin 2 + n\ 2 
Using (|8| we calculate the density of the corresponding packing, 

a+1 7T (a + l)(nl + riin 2 + n 2 ) 



6 



7T 

4 



n|+nin 2 +n| 



(a + l)(nf + n x n 2 + nf) 
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Since the packings coming from strip tilings are not a triangular lattice pack- 
ing 5 < 5a, and thus 



c > (a + l){n\ + nin 2 + > "a + ni^ 2 + nj, (11) 

and 

n\ + nin 2 + n\ ^ 
c 

Then we can improve the inequality (11) to 

n / 9 9s 7i? + nino + n 2 . 
c > (a + l)(n? + 7ii7i 2 + nl) + ^ - 

o / \/Q *3 \ *3 Q 

c > c(o + l)(n 1 + riin 2 + n 2 ) + «i + niri2 + n 2 
c 2 > [c(a + l) + l}(n 2 1 + n 1 n 2 + n 2 2 ). 

Then going back to the total number of packing elements n = c(a +1), we 

get 

n c(a + l) (a + l)(n? + riin 2 + n|) 



Ti+1 c(a + 1) + 1 



> 



Thus 



n . . (a + l)(n? + n x n 2 + n%) . 
-oa > oa = o. 



n + 1 c 

□ 

Note that it is not used that such triangle-rhombus strip packings are 
strictly jammed, but they are. It is easy to find a self-stress on the (strut) 
framework. For the edges that are in the direction of g& choose a positive 
stress ui, and for the direction g\ choose u 2 and for the direction g choose 
u 3 such that u>ig& + oj 2 g\ + g) = 0. Then any constant positive stress 
on the horizontal members will do. 

When a > 0, the connected component of triangles is infinitesimally rigid 
all by itself, since it is connected in the torus. The case when a = was 



essentially handled in Subsection 4.1 and Equation 



So Theorem 10 verifies the density gap conjecture Conecture[3]in the case 
of strips of rhombi and triangles. 

In order believe that these calculations correspond to actual packings on 
a triangular torus, we present a specific example for all values of a = 0, 1, ... . 
For any non-negative integer a, let b = 4a + 5, c = (a + 1)7 + 1 = 7a + 8, 
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Tlx = 2, n 2 = 1, Us = 1, and n 4 = 1. Then 



nin 4 - n 2 n 3 
n\ + riin 2 + n 2 2 



9 



1 

7 

7a + 8 



7 



(2 + Sa 1 ; 



1 
7 



8>/3. 
7~T 



One can verify that Equation ( 10 ) holds for any a > 0. Thus this corresponds 
to collectively jammed strip packings. Figure [TT] corresponds the case when 
a = 0. 



5 Relation to prevoius work 

One of the starting points for this work was from an example of a packing 



provided by Ruggero Gabbrielli who suggested that the packing in Figure [13 
was at least collectively jammed. At first it was difficult to determine if, for 
this this packing, it was true, but it was shown by the first author at a talk 
at the Fields Institute |3] that it was not collectively jammed. At the next 
talk by Walter Whiteley [20] it was shown that by looking at parallel classes 
of edges of a graph one can create a flex, which is a basic idea of Theorem 
|5} It was also shown in [IT] . Section 4.9 and Figures 21-22 how the rotating 
motion described work, and it was shown that there is such a flex when there 
are two classes of edges as there is here. 

Triangle- square tilings have been of some interest in the physics literature 
as in [TU [15J Uni EE]. Indeed, the idea of associating to a triangle-square 
tiling, a graph, where every vertex is of degree six, and identifying it with 
a triangle lattice, is in [H]. In [13] and [13] they estimated the ratio of 
squares to triangles in a "random" square-and-triangle tiling. Square lattice 
tilings were quite popular. Suppose that that one side of a period square of 
square lattice has a polygonalpath involving the vectors 1, g 2 , g, g~ x , where 
g = v / 3/2 + i/2 as in Section |2j So the sum of the edge vectors that go from 
one corner to the next horizontal corner of the square will be 

s = a + b(g + g^ 1 ) = a + b\/3, 

where a and b are integers, since 1, g 2 , g, g~ l is another basis for the lattice 
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Figure 13: This is a packing in a triangular torus with a 12-gon face that, 
for a generic configuration of the packing graph is rigid, but it is not rigid as 
a packing graph since it is a subset of the packing in Figure |4j 



determined by the edge vectors. But the area of the square torus is 

s 2 = a 2 + 3b 2 + 2abV3 = / 4 + ^/ 3 , 

where f% is the number of triangles and / 4 is the number of squares. So 

/ 3 = 8ab and / 4 = a 2 + 3b 2 . 

Thus 

ft a 2 + 3i> 2 1 m 3 ^ 

since the minimum of a; + 3/a; for x > is when x = v3 as any calculus 
student knows. 

A similar analysis for the triangular torus shows that y/3/<i is a maximum 
for the ratio fi/fz- 

The ratio V3/4 for f±/ f$ was claimed as the maximum entropy for triangle- 
square tilings with twelve-fold symmetry in |13j . 
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6 Future work 



It would be helpful to understand the case when the graph of the packing 
consists of triangles and rhombi more generally than in Theorem [6] There is 
the question of what kind of packing graphs can occur. 

More ambitiously it would be very useful to understand what kind of 
collectively jammed packing graphs are that have faces of degree larger than 
four. One idea is to complete the graph of the packing to a triangulation, 
where the degree of each vertex is six. Then there are ways to compare 
the density of such packings to the density <5a and ^x^a using some simple 
techniques from the theory of the global rigidity of tensegrities. 
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